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ABSTRACT 

The results of three-dimensional numerical simulations of mass transfer in binaries 
with spiral shocks are presented. It is shown that the mass transfer rate variation 
disturbs the equilibrium state of accretion disk and results in formation of the blob. 
This dense formation moves on the disk with variable velocity. The results of our 
simulations show that the blob lives long enough and retains its main characteristics 
for the time of order of tens orbital periods. Analysis of the results shows that the 
action of dissipation is negligible on these timescales and the blob smearing out due 
to differential rotation of the disk is stopped by its interaction with spiral shocks. 

Light curves of some CVs show periodic or quasi-periodic photometric modula- 
tions with typical width ~ r !l^0 r !2 (i.e. with period ~ 0.1-^0.2P or (,). We suggest to 
consider the formation of the blob in binaries with spiral shocks as a possible reason 
for these modulations. The efficiency of this mechanism is confirmed by good quali- 
tative (and, in part, quantitative) agreement between the results of our simulations 
and observations of binaries IP Peg and EX Dra. 



1 Introduction 

Analysis of light curves in CVs with high or- 
bit inclination and especially eclipsing bina- 
ries permits to get information on the flow 
structure between the components of the sys- 
tem. The shape of light curves in these bina- 
ries is very sophisticated and even uneclipsed 
parts of light curve show significant bright- 
ness .oscillations (see, i-e.g-, Hack & La Dous 
1993G3; Warner 19950). A part of these os- 
cillations, the so called flickering has aperiodic 
nature and is characterized by small magni- 
tude and short timescale. On the other hand, 
besides sporadic short-period oscillations light 
curves in CVs demonstrate periodic or quasi- 
periodic photometric modulations with typi- 



cal period ~ 0.1 -r- 0.2P or b. Various models 
were suggested to explain these modulations 
of light curves: from intrinsic disk instabil- 
ities up to modulations of the rate of mass 
transfa from the mass-donating star (Warner 
1995EI). Based on the results of 3D gasdy- 
namical simulations we concluded in (Bisikalo 
et al. 2001 0, hereinafter Paper I) that if spiral 
shocks are presented in the accretion disk then 
any disturbance of the disk would result in the 
appearance a new dense formation in the disk 
- a blob, the latter moving through the disk 
with variable velocity. A significant density 
contrast between the blob and the disk (up to 
1.5 times) as well as the variable velocity of 
the blob (but with constant period of revolu- 
tion ~ 0.1 ^-0.2P or b) permit us to consider the 
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appearance of this formation as a possible rea- 
son for observable quasi-periodic oscillations 
on the light curves in CVs. If the cmissivity of 
the blob is larger than that for the accretion 
disk (e.g., for the case of recombination radi- 
ation which is proportional to the square of 
density) the oscillations on the light curve can 
be caused by the occultation the blob by op- 
tically thick disk. And if the blob emissivity 
is less than that for the disk the oscillations 
can be caused by occultation of the disk by 
optically thick blob. 

The main aim of this work is more detailed 
(as compared to Paper I) investigation of the 
blob formation as well as a determination of 
the dependence of its features on the param- 
eters of the model. We also analyze observa- 
tions of eclipsed binaries with spiral shocks (IP 
Peg and EX Dra) in the context of applicabil- 
ity of proposed mechanism for explanation of 
oscillations on light curves in these binaries. 

2 Results of numerical 

simulations in CVs with 
variable mass transfer 
rate 

In Paper I we have presented the results of 3D 
gasdynamical simulations of flow structure in 
semidetached binaries after the mass transfer 
termination. The investigation of the resid- 
ual accretion disk shows that as early as at 
the time t = <o + 0.2P oro (where to is time 
when mass transfer has been ceased) the flow 
structure is changed significantly. The stream 
from L\ doesn't dominate anymore, and the 
shape of accretion disk changes from quasi- 
elliptical to circular. The second arm of tidally 
induced spiral shock is formed while earlier 
(before the termination of mass transfer) it 
was suppressed by the stream from L±. The 
results of numerical simulations presented in 
Paper I show that a dense blob is formed in 
the residual disk, the velocity of the blob mo- 
tion through the disk being variable. The blob 
is a dense formation which in steady-state so- 
lution coincides with a post-shock zone of the 
spiral shock. After mass transfer termination 
when the nourishment of the disk from the 
stream is ceased and equilibria of the disk is 
disturbed the blob comes off from the front of 



spiral shock and begins to move. 

Two factors can force the blob to smear 
out as it moves through the accretion disk 
- dissipation and differential rotation of the 
disk. Dissipation smearing occurs in viscous 
timcscalc which can be evaluated as tens - 
hundreds orbital periods for a ~ 0.01 -j- 0.1 
which is typical for CVs.[] We would like to 
stress that the viscosity timescale of smear- 
ing out of the blob gives the lower estimation 
of its lifetime, since this time is so large that 
the recurrent mass transfer rate variation and 
the subsequent disturbance of the disk is very 
probable. It means that the new disturbance 
of the disk due to mass transfer rate variation 
can occur earlier than the dissipative smearing 
out of the blob. 

Smearing out (or mixing) of the blob due 
to differential rotation of the disk is virtu- 
ally more important since this process occurs 
in dynamical timescale and one may suppose 
that the blob disappears during few orbital 
periods. Nevertheless, 3D gasdynamical cal- 
culations of Paper I show that existence of 
a spiral shock prevents dynamical mixing of 
the blob. Passing through the front of spi- 
ral shock the blob becomes denser and more 
compact, this effect being presented up to the 
moment when spiral shock disappears. The 
density contrast between the blob and the disk 
begins to decrease after the disappearing of the 
spiral shock as well. 

We understand that the situation of com- 
plete termination of mass transfer from donor- 
star is rather rare, so we should consider the 
possibility of the blob formation when the vari- 
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where Rd — the radius of accretion disk; v - kinematic 
viscosity; c - sound speed; a - Shakura-Sunyaev pa- 
rameter; H — height scale of the disk; /3 = H/R d = 
c/Vk', &k and Vk — angular and linear rotational ve- 
locities of Keplerian disk. Transform the last expres- 
sion to obtain 

tv = n orb = ^T+g (Rd\ 3/2 

Port 2naj3 2 n K 2naj3 2 \ A J 

where P orb - orbital period; q - mass ratio (the ratio 
of donor-star mass to accretor mass); A - the separa- 
tion of binary. In our calculations the typical values of 
parameters are: /3 ~ 0.1 -f 0.15 and Rd/A ~ 0.3 4- 0.4, 
so for a ~ 0.1 we have t„ ~ 25P or (, while for a ~ 0.01 
we have r„ ~ 250P o rt- 111 Paper I we have obtained 
the lifetimes of accretion disk for various values of a 
which give a little bit less but comparable values of t„ . 
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ation of the mass transfer rate is relatively 
small. In this paper we present the results of 
investigation of the flow structure in semide- 
tached binary after the mass transfer becomes 
weaker twice and 10 times. 

Note, that in the case of complete termi- 
nation of mass transfer our calculations have 
been stopped after the disk vanishes, while 
in case of the reduced rate of mass transfer 
we have observed the flow structure went out 
to a new quasi-steady-state solution. In gen- 
eral, to reach an ultimate steady-state solu- 
tion we should conduct our simulations over 
time interval comparable to viscous timescale 
of the disk. Nevertheless, we have restricted 
this time to few tens of orbital periods bear- 
ing in mind the typical timescale of outburst 
activity in CVs resulting in the formation of 
spiral shocks. 

The full description of the 3D gasdynami- 
cal model can be found in Paper I. Here we 
pay attention only to the main features of the 
model. To describe the gas flow in the bi- 
nary we used the 3D system of Euler equa- 
tions for Cartesian coordinate system. To 
close the system of equations, we used the 
equation of state of ideal gas with adiabatic in- 
dex 7 ~ 1, that corresponds to the case close 
to the isothermal one. Other dimensionless 
parameters which determine the solution were 
adopted as follows: mass ratio is q = 1, the 
ratio of sound speed in L± to orbital veloc- 
ity is e = c(Li)/Att. The gas flow was simu- 
lated over a parallelepipedon [ 1 / 2 A . . . 3 / 2 A] x 
[- 1 / 2 A . . . 1 / 2 A] x [0 . . . 1 /iA] without a sphere 
with a radius of / ioo^4 representing the accre- 
tor. The calculations were conducted on the 
grid 61x61x17, in terms of the a-disk the nu- 
merical viscosity approximately corresponded 
to a ~ 0.05. To obtain the numerical solu- 
tion of the system of equations we used the 
Roe-Osher TVD scheme of a high approxima- 
tion order (Roe 19860 ; Chakravarthy & Os- 
her 19850). We used quasi-steady-state solu- 
tion for constant rate of mass transfer as initial 
condition. At the moment of time t = to we 
decreased the density of the injected matter 
twice (run "A") and 10 times (run "B"). 

The Figure 1 presents time variation of the 
mean density of matter passing through a 
semi-plane XZ (Y = 0, X > A) slicing the 
disk. Along with a general density fall we 
can see here periodic oscillations called forth 
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Figure 1 : Time variation of the mean density of mat- 
ter passing through a semi-plane XZ (Y = 0, X > A) 
slicing the disk for run "A" (upper curve) and "B" 
(lower curve). 



by passing the blob. The analysis of curves 
in Fig. 1 shows that after changing (decreas- 
ing) the mass transfer rate a dense formation 
- blob - arises. The period of blob revolu- 
tion is ~ 0.15P or f, for run "A" and ~ 0.18P or b 
for run "B" . After recurrent reaching of quasi- 
steady-state solution (i.e. for t > to + 3P or b for 
run "A" and t > t + 13P or6 for run "B") the 
period of the blob revolution increases and be- 
comes ~ 0.16P orb for run "A" and ~ 0.20P orb 
for run "B". The new quasi-steady-state so- 
lution is characterized by long-period oscilla- 
tions (Fig. 1) with period ~ 5.1P orb for run 
"A" and - 2.5P orb for run "B". The Fig- 
ure 2 presents the Fourier analysis of curves 
from Fig. 1. It is seen that for run "A" domi- 
nating harmonics are 0.155P or f, and 0.165P orb 
(Fig. 2a). The beating of these two gives long- 
period oscillations with period ~ 5.1P orb . For 
run "B" dominating harmonics are 0.18P orb 
and 0.21P orb (Fig. 2b), and their beating 
gives long-period oscillations with the period 
~ 2.5P orb . 

Let us consider how does the disk struc- 
ture change after the rate of mass transfer 
decreases. The Figure 3 presents the re- 
sults of run "A" and shows the density dis- 
tributions over the equatorial plane for 6 mo- 
ments of time t + 0A4:P orb , 0.47P or6 , 0.49P orb , 
Q.52P orb , 0.54P orb , 0.57P or6 , covering the full 
period of the blob revolution. The maxi- 
mal density in Fig. 3 corresponds to p ~ 
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Figure 2: Strength of Fourier harmonics versus pe- 
riod for run "A" (upper panel) and "B" (lower panel). 

0.035p(Li). These results are relevant to the 
transient stage when dense stream from L\ 
vanishes at the moment of time t = to & n d 
the flow structure is determined by the rela- 
tively dense disk while the matter of rarefied 
stream (with density twice less than the ini- 
tial density) plays a minor role. This stage 
begins after the last portions of matter of the 
old stream come to the disk, i.e. at t = 
to + 0.2P or b, and finishes when a new quasi- 
steady-state solution will be reached. The 
main features of this stage are: the dense resid- 
ual disk dominates; its shape changes from 
quasi-elliptical to near circular; the second 
arm of spiral shock forms. 

The Figure 4 also presents the results of 
run "A" and shows the density distributions 
over the equatorial plane for other 6 moments 
of time t + 17.64P orf) , 17.67P or6 , 17.70P or6 , 
17.73P or6 , 17.76P orb , 17.79P or6 , covering the 
full period of the blob revolution for the stage 
when the new quasi-steady-state solution is 
established. The maximal density in Fig. 3 
corresponds to p ~ 0.02 j o(ii). The analy- 
sis of these drawings shows that obtained flow 
structure is similar to the solution for the case 
of constant rate ofjnass transfer (see, e.g., 
Bisikalo et al. 1998^, 20000, Paper I). It is 



seen that the stream from L\ dominates. It is 
also seen that the uniform morphology of the 
system stream-disk results in quasi-elliptical 
shape of the disk and the absence of 'hot spot' 
in zone of stream-disk interaction (X ~ 0.8A, 
Y ~ —0.1 A). At the same time the interaction 
of the gas of circumbinary envelope with the 
stream results in the formation of an extended 
shock wave located along the stream edge ('hot 
line'). The Figure 4 also manifests the forma- 
tion of only one arm of tidally induced spiral 
shock located in I and II quadrants of coor- 
dinate plane. The flow structure in the place 
where the second arm should be is defined by 
the stream from L\ which appears to suppress 
the formation of another arm of the tidally in- 
duced spiral shock. 

Nevertheless, opposite to the case of con- 
stant rate of mass transfer the blob remains 
to live and gives the periodic changes of disk 
structure. The Figure 4 shows that the blob 
tends to spread uniformly but passing through 
the arm of spiral shock it is retarded, this leads 
to the formation/sustaining of compact blob. 
This also leads to periodic changes of disk 
structure in a new quasi-steady-state solu- 
tion. The Fourier analysis shows (see Fig. 2a) 
that for run "A" dominating harmonics are 
0.155P or & and 0.165P or &. The beating of these 
two gives long-period oscillations with period 
~ 5.1P rb- First (short-period) harmonic ap- 
pears to be due to interaction of the blob 
with the arm of spiral shock and the second 
harmonic is due to interaction of peripherical 
parts of the blob with the stream of matter 
from Li, the period of latter harmonic being 
larger. The larger period is result of differen- 
tial rotation of the disk - peripherical parts 
of the blob locate at larger radii and rotate 
slower. After interaction with the stream this 
lag decreases and the blob becomes compact. 

Comparative analysis of Figs 1-4 shows that 
the mean period of the blob rotation is the 
same both for transient stage when two arms 
of spiral shock exist and for new quasi-steady- 
state solution with one arm of spiral shock. 
At the same time the density contrast between 
the disk and the blob changes significantly due 
to long-period modulations (Fig. 5). It is seen 
that on the transient stage when the dense disk 
dominates the contrast is ~ 1.5. Later, after 
reaching the recurrent quasi-steady-state solu- 
tion with two dominating harmonics and long- 
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Figure 3: Density distribution over the equatorial plane for run "A". Results are presented for the moments 
of time to + 0.44P ori) , 0A7P orb , 0A9P orb , 0.52P orh , 0.54P ori) , 0.57P or6 , covering the full period of the blob 
revolution for the transient stage. The maximal density corresponds to p ~ 0.035p(Li). 
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Figure 4: The same as in Fig. 3 for other six moments of time: tQ-\-\l .64P or b, 17.67P or ^, 17.70-P or ^, 17.73-P or ^, 
17.76P or (,, 17.79P or (,, covering the full period of the blob revolution for the stage when the new quasi-stcady- 
state solution is established. The maximal density corresponds to p ~ 0.02p(Li). 
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Figure 5: The density contrast between the disk and 
the blob versus time. 



period modulations the density contrast varies 
from ~ 1.1 to - 1.8. 

Resuming the results of run "A", we can 
conclude that the solution at the recurrent 
quasi-steady-state stage is defined by the 
stream from L\ and the presence of the blob 
manifests in periodic variations of the disk 
structure. 

Let us consider now the results of run "B" 
when the mass transfer rate was decreased 
10 times and the gas of the stream has sig- 
nificantly less density. Similar to run "A" 
on the transient stage the dense disk and 
the blob dominate. But opposite to run 
"A" the blob continues to dominate even at 
the recurrent quasi-steady-state stage. Oppo- 
site to run "A" the blob changes the general 
flow structure, in particular, the blob 'blows 
off' the stream when reaching it. This is 
seen from Fig. 6 where the density distribu- 
tions over the equatorial plane for moments 
of time t + 22.45P orb , 22.48P orb , 22.51P Q r&, 
22.54P or6 , 22.57P or6 , 22.61P orf) , covering the 
full period of the blob revolution for the new 
quasi-steady-state stage of run "B" . The max- 
imal density in Fig. 6 corresponds to p ~ 
0.0035p(Pi). The small value of density con- 
firms the significant influence of the blob on 
the flow structure on all stages of run "B". 
The Fourier analysis of results of the run "B" 
shows (see Fig. 2b) that dominating harmon- 
ics are 0.18P or b and 0.21P or (,. The beating of 
these two gives long-period oscillations with 
period ~ 2.5P or b- Similar to run "A" these 



harmonics appears to be due to interaction of 
the blob with the arm of spiral shock and with 
the stream of matter from L±. 

The analysis of these results as well as the 
results of Paper I shows that for binaries with 
spiral shocks small disturbances resulting from 
variations of mass transfer rate lead to forma- 
tion of the blob. Moreover, on the dynamical 
timescale the blob is supported by the pres- 
ence of spiral shocks and this mechanism is 
independent on the variations of mass trans- 
fer rate. Thus, we can conjecture that not only 
mass transfer rate variation but every distur- 
bance of disk structure will transform into the 
blob. The lifetime of the blob is determined by 
the lifetime of spiral shock which prevents the 
blob spreading in dynamical timescale, and 
dissipation which governs the blob spreading 
in viscous timescale. 

3 The blob parameters and 
the variability of light 
curves of CVs with spiral 
shocks 

Spiral shocks or, more exactly, the dense post- 
shock regions attached to the arms of the 
spiral shock were often considered as proba- 
ble source of variability of astrophysical ob- 
jects (see, e.g., Chakrabarti & Wiita 1993a0, 
1993btf). Recently we argued in Paper I that 
the dense blob can detach from the front of 
spiral shock and further travel through the 
disk. In Paper I as well as in the present work 
the blob is formed as a result of disturbance 
of the disk structure due to variation of mass 
transfer rate. Once forming blob tends to dis- 
tribute uniformly along the disk due to differ- 
ential rotation but retards after passing trough 
spiral shock so the compactness of the blob 
retains/resumes. Later the increasing den- 
sity contrast between the blob and the disk as 
well as growths of pressure gradient will force 
the standstillcd blob to move. When reaching 
again the front of spiral shock the process of 
formation/rctainment of the blob is repeated. 
Results of 3D gasdynamical calculations show 
that the blob travels through the disk during 
tens of orbital periods, its main characteris- 
tics being conserved. The blob ought spread 
under the action of dissipation but the vis- 
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Fi gurG 6: Density distribution over the equatorial plane for run "B". Results are presented for the moments 
of time t + 22.45P ori) , 22.48P ori) , 22.51P or6 , 22.54P or6 , 22.57P or6 , 22.61P or6 , covering the full period of 
the blob revolution for the new quasi-steady-state stage of run "B". The maximal density corresponds to 
p ~ 0.0035p(Li). 
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cous timescale is much greater than the time 
of our calculations so we can neglect the vis- 
cous spreading of the blob. We should like to 
stress that the we have chosen the duration 
of our calculations (up to ~ 25P or b) not ar- 
bitrary but in agreement with typical length 
of the outburst decline (tens of P or b)- It is 
on the stage of outburst decline the decreas- 
ing of mass transfer rate is possible, and, be- 
sides, on this stage the spiral shocks are ob- 
servable (or, more exactly, the most reasonable 
assumptions on their existence can be made 
here) . 

To determine the typical features of the blob 
in CVs it is necessary to conduct a set of 
calculations covering all the range of param- 
eters which are relevant to CVs. Neverthe- 
less, based on the results of Paper I as well as 
on the results of this work we can conclude 
that generally the period of the blob revo- 
lution depends on the value of viscosity: it 
equals to 0.15 0.20P or b for a ~ 0.05, and to 
0.10 0.15P or6 for a ~ 0.01. This means that 
for typical accretion disk with a ~ 0.01 -j- 0.1 
the period of the blob revolution should be 
in range 0.10 -j- 0.25P O r&- Another important 
property of the blob is the density contrast 
w.r.t. the rest of the disk since this parame- 
ter determines the observational evidences of 
the blob. As it was shown above the value of 
density contrast is independent on the rate of 
mass transfer and changes periodically around 
the mean value of ~ 1.5. The nature of mech- 
anism of formation/sustaining of the blob per- 
mits us to believe that this value is typical for 
different CVs. 

Let us consider the virtual observation ev- 
idences of the blob and compare it with the 
real observations. The existence of the blob 
is maintained by spiral shocks, so we have 
chosen the observation of IP Peg which was 
gave the first evidences of the spiral shocks 
(Steeghs, Harlaftis & Home 1997E3 ; Harlaftis 
et al. 19990) and EX Dra where the exis- 
tence of spiral shocks is also, probable (Bap- 
tista Catalan 19990, 2000i; Joergens et,al. 
2000E3; Joergens, Spruit & Rutten 2000113). 

Since the blob is preserved from the spread- 
ing by the existence of spiral shock let us con- 
sider the link between the light curve vari- 
ations and the appearance/disappearance of 
spiral shocks. The standard treating of IP Peg 
and EX Dra suggests the extension of the ac- 



cretion disk and formation of spiral shocks in 
it at some stage. After the outburst ends and 
the quiescence state establishes the radius de- 
creases so tidal interaction can't induce spi- 
ral shocks anymore. Thus, the light curves 
should manifest quasi-periodic variations only 
during outburst and, possible, some time af- 
ter it.Q Observations show variations of un- 
eclipsed parts of the light curves for both sys- 
tems during outburst and immediately after it. 
The width of these variations is ~ 0^1-^0^2 for 
EX Dra (Baptista, Catalan & Costa 2000H) 
and for IP Peg CMarsh & Home 1990O; 
Webb et al. 199913), the amplitude of vari- 
ations sufficiently decreasing when the quies- 
cence state establishes. These observations 
confirm the proposed model on the qualitative 
level: it is seen that variations of uneclipsed 
parts of light curves correlate with the pres- 
ence of spiral shocks, and the width of varia- 
tions corresponds to the theoretical one. 

As another test for the suggested model let 
us consider the quantitative evidences of the 
blob existence. The analysis of light curves 
shows that the disk deposits the major fraction 
of luminosity during outburst (up to_?/ 6 for IP 
Peg, see, e.g., Khruzina et al. 2001 liJ), so the 
amplitude of light curve variations should be 
approximately equal to the density contrast. 
Our calculations give the mean value of den- 
sity contrast as 1.5 so taking into account the 
relative luminosity of the disk we can expect 
the amplitude of variations on light curves as 
~ 40%. Observations show the amplitude of 
variations of the flux on uneclipsed parts of 
light curves as ~ 30% for IP Peg (see, e.g., 
Steeghs et al. 1996^3; Morales- Rueda, Marsh 
& Billington 20000) and - 15%, for EX Dra 
(Baptista, Catalan & Costa 2000B). Unfortu- 
nately, the analysis of light curve variations in 
quiescence is inconvenient due to both decreas- 
ing of density contrast and relative luminosity 
of the disk ( 1 / 5r ipr IP Peg, see, e.g., Khruz- 
ina et al. 2001 EJ), so we restrict ourself by 
consideration of quantitative estimations for 
outburst only. The comparison of these esti- 
mations confirms the proposed model on the 
quantitative level as well. 

As an additional test of the model applica- 

^ote that if the viscosity in the disk is small then 
the dissipative spreading is negligible and the blob can 
exist even when the spiral shocks vanish. Correspond- 
ingly, the light curve variations can be observable after 
that while with decreasing magnitude. 
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bility we would like to stress its relevance to of President of Russia (99-15-96022, 00-15- 
the explanation of variation of emission lines 96722). Authors wish to thank A.M.Fridman 
profiles in CVs. For example, the observations and O.V.Khoruzhii for useful discussions, 
of H a profiles in the disk_af binary SS Cyg 
(Martmez-Pais et al. 1994E9) show its essen- 
tially asymmetrical nature, the distance be- 
tween peaks being significantly variable. The 
presence of spiral shocks permits to explain 
the asymmetry ppf line profiles (Chakrabarti 
& Wiita 1993atP), but the spiral shocks by 
itself can't give the mechanism resulting in 
variable distance between peaks. The mecha- 
nism suggested in this work do explain both 
the asymmetry of line profiles (in the same 
manner as above) and the variability of the 
distance between peaks (by the variability of 
the projection of blob velocity). Note also 
that the sporadic low-magnitude variations on 
the light curve (flickering) can also be a con- 
sequence of the presence of the blob in the 
disk. Disturbances of the disk structure due 
to periodic travelling of the blob as well as due 
to stream/blob interaction can result in low- 
amplitude oscillations manifesting as sporadic 
variations on the light curve. 

4 Conclusions 

Light curves of CVs show periodic or quasi- 
periodic photometric modulations with typi- 
cal width - 0^1 ~ P 2. The 3D gasdynam- 
ical simulations of flow structure in binaries 
disturbed by the mass transfer rate variation 
prove the appearing a new dense formation - 
blob, the latter moving through the disk with 
variable velocity. Long lifetime of the blob, 
significant density contrast between the blob 
and the rest of accretion disk (the mean value 
is ~ 1.5), and characteristic period of its rev- 
olution ~ 0.1 -T- 0.2P or b permit us to consider 
this formation as a possible reason for observ- 
able quasi-periodic variations of light curves 
of CVs. Note that if suggested mechanism is 
valid then the variations of light curves of CVs 
give an additional proof for the spiral shocks 
existence. 
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